NOBO_GoalSetting

NOBO interactive Map: GCT & Trend Data

This document creates a leaflet map showing goals collected as part of the Northern Bobwhite, Grasslands and Savannas National Partnership along side data collected showing NOBO trends collected by James Martin

Loading Libraries

library(tidyverse)
library(sf)
library(leaflet)

## Warning: package 'leaflet' was built under R version 4.0.3

library(htmltools)

## Warning: package 'htmltools' was built under R version 4.0.5

library(htmlwidgets)
library(raster)
library(gstat)

## Warning: package 'gstat' was built under R version 4.0.3

library(spatial)
library(dplyr)
library(jsonlite)

## Warning: package 'jsonlite' was built under R version 4.0.5

library(ggplot2)
library(hrbrthemes)

## Warning: package 'hrbrthemes' was built under R version 4.0.5

library(ggthemes)

options(scipen=999)

Loading Variable Sets

##Loading Boundaries
State_Boundaries <- st_read("cb_2018_us_state_500k.kml", quiet = TRUE)
State_Boundaries_GJSON <- st_read("US_State_Boundaries.json", quiet = TRUE)
State_Boundaries_Zip <- st_read("C:/Users/sageg/Desktop/newRrepo/USDA_WLFW_NOBO/StateBoundariesZIP", quiet = T)



##Loading GCT Data
GCT_data <- read_csv("NOBODATA_ForLeaflet_Final.csv")

##Loading RDS
bird_data <- readRDS("NOBO_route_level_trends.rds")

###Adjusting NOBO Trend Datasets

separated_coord_bird_data <- bird_data %>%
    mutate(lat= unlist(map(bird_data$geometry,1)),
           long= unlist(map(bird_data$geometry,2)))

###Adding Geometries to GCT Data

##Joining Boundaries and GCT Data
GCT_and_Geographies <- st_as_sf(left_join(GCT_data, State_Boundaries_Zip, by =c("STATE" ="NAME")))

Creating Label for GCT Data

## Creating Label
GCT_and_Geographies$label <- 
  paste("<b>", "<big>", GCT_and_Geographies$STATE,"</b>", "</big>",
       "<br>", GCT_and_Geographies$REGION,
        "<br>", "Values shown below reflect the combined value of", "<br>", "Original Program Goals and Framework Goals",
       "<br>",
        "<br>",
        "<b> Top 3 Core Conservation Practices: </b>",
       "<br>",
        GCT_and_Geographies$TOP3CORE,
        "<br>",
       "<br>",
        "<b>", "Financial Assistance:", "</b>","<br>",
        "$ ", prettyNum(GCT_and_Geographies$REQFIN_FRAME,big.mark=","),
       "<br>",
        "<br>",  "<b>", "Total CP Coverage, Acres:", "</b>",
        "<br>", "Core:", prettyNum(GCT_and_Geographies$ACRE_FRAME_CORE, big.mark = ","),
        "<br>", "Supplemental:", prettyNum(GCT_and_Geographies$ACRE_FRAME_SUPP, big.mark = ","),
        "<br>", "Core and Supp Combined:", prettyNum(GCT_and_Geographies$ACRE_FRAME_CAS, big.mark = ","),
       "<br>",
        "<br>", "<b>", "Total CP Coverage, Feet:", "</b>",
        "<br>", "Core:", prettyNum(GCT_and_Geographies$FT_FRAME_CORE, big.mark = ","),
        "<br>", "Supplemental:", prettyNum(GCT_and_Geographies$FT_FRAME_SUPP, big.mark = ","),
        "<br>", "Core and Supp Combined:", GCT_and_Geographies$FT_FRAME_CAS,
       "<br>",
        "<br>", "<b>", "Total CP Coverage, Number of X:", "</b>",
        "<br>",GCT_and_Geographies$X_FRAME,
       "<br>",
        "<br>", "<b>", "Number of Written Plans:", "</b>",
        "<br>",GCT_and_Geographies$WRITTEN_FRAME,
       "<br>",
        "<br>", "<b>", "Number of Applied Plans:", "</b>",
        "<br>",GCT_and_Geographies$APPLIED_FRAME) %>% 
  lapply(htmltools::HTML)

###Setting up Leaflet with only GCT Data

nobo_summary_leaflet <- leaflet(GCT_and_Geographies)%>%
  addProviderTiles(providers$Stamen.Watercolor) %>%
  addPolygons(highlightOptions = highlightOptions(fillOpacity = 1, fillColor="cornsilk"),
              weight=1, 
              fillColor= "darkseagreen",
              color = "black",
              fillOpacity=.5,
              label = ~GCT_and_Geographies$label, 
              labelOptions = labelOptions(
    style = list("font-weight" = "normal", padding = "3px 8px"),
    textsize = "11px",
    direction = "auto")) %>%
  
    addControl("Northern Bobwhite, Grasslands, and Savannas National Partnership 2022-2026: Framework Goal Summary", position = "topright")

## Warning: sf layer has inconsistent datum (+proj=longlat +datum=NAD83 +no_defs).
## Need '+proj=longlat +datum=WGS84'

nobo_summary_leaflet 

saveWidget(nobo_summary_leaflet, file="NOBO_SUMMARY_LEAFLET.html")

leaflet with both GCT and Trend Data

## Setting up LEaflet 
nobo_summary_leaflet <- leaflet()%>%
  addProviderTiles(providers$Stamen.Watercolor) %>%
  addPolygons(data=GCT_and_Geographies, 
              highlightOptions = highlightOptions(fillOpacity = 1, fillColor="cornsilk"),
              weight=1, 
              fillColor= "darkseagreen",
              color = "black",
              fillOpacity=.5,
              label = ~GCT_and_Geographies$label, 
              labelOptions = labelOptions(
    style = list("font-weight" = "normal", padding = "3px 8px"),
    textsize = "11px",
    direction = "auto")) %>%
  
  ## Add lci-data
  addCircles(data=separated_coord_bird_data,
            lng = ~long, 
            lat = ~lat)%>%
  
    addControl("Northern Bobwhite, Grasslands, and Savannas National Partnership 2022-2026: Framework Goal Summary", position = "topright")

## Warning: sf layer has inconsistent datum (+proj=longlat +datum=NAD83 +no_defs).
## Need '+proj=longlat +datum=WGS84'

nobo_summary_leaflet 

#saveWidget(nobo_summary_leaflet, file="NOBO_SUMMARY_LEAFLET.html")

Showing BirdData as GGPLOT and Leaflet

## ggplot
ggplot(bird_data, aes(color=lci_trend)) +
  geom_sf()

## leaflet

leaflet(separated_coord_bird_data) %>%
addProviderTiles("CartoDB.Positron") %>%
  addCircles(lng = ~long, 
            lat = ~lat, 
            color ="red")